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1.  Abstract 

Bypass  transition  of  a boundary  layer  takes  place  at  high  main  stream  tur- 
bulence levels.  It  is  governed  by  the  intrusion  of  non-linear  disturbances  into 
the  viscous  sublayer.  In  this  paper  a numerical  study  is  presented  in  which  a 
spatially  developing  laminar  boundary  layer  is  subdued  to  large  scale  (non- 
linear) disturbances.  The  (fully  compressible)  Navier-Stokes  equations  are 
solved  using  direct  numerical  simulation  with  a higher  order  finite  difference 
method  on  a collocated  grid.  The  results  show  that  the  initially  smooth  dis- 
turbance rapidly  changes  into  an  arrow-head  structure  with  the  characterics 
of  a turbulent  spot. 

2.  Introduction 

In  the  design  of  turbomachinery  the  prediction  of  the  boundary  layer  layer 
heat  transfer  is  of  mayor  importance,  to  since  the  hot  gases  from  the  com- 
bustion chamber  force  an  enormous  heat  flux  to  the  turbine  blades.  The 
heat  transfer  in  a boundary  layer  increases  significantly  when  laminar  to 
turbulent  transition  occurs.  When  transition  starts  it  is  assumed  that  so 
called  turbulent  spots  are  initiated  in  the  laminar  boundary  layer. 

Turbulent  spots  have  been  first  observed  by  (Emmons,  1951)  in  a water 
channel.  Most  transition  models  which  are  used  nowadays  are  based  on 
these  turbulent  spots.  It  is  supposed  that  a turbulent  spot  is  a local  area 
of  turbulence,  that  convects  in  streamwise  direction  with  a mean  velocity 
which  is  less  than  the  main  flow  velocity.  The  leading  edge  travels  faster 
than  the  trailing  edge  and  so  the  turbulent  spot  grows  lengthwise  while 
moving  in  downstream  direction.  Furthermore,  the  spots  appear  to  remain 
at  the  same  shape  so  they  also  grow  laterally.  At  a certain  point  the  spots 
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start  to  merge  untill  the  flow  is  completely  turbulent.  This  is  said  to  be  the 
end  of  transition. 

In  general  a distinction  between  natural  and  bypass  transition  must  be 
made.  At  low  disturbance  levels  in  the  main  flow  Tollmien-Schlichting  waves 
are  initiated  at  a certain  Reynolds  number.  These  waves  grow  in  amplitude 
and  three  dimensional  waves  start  to  develop.  Further  downstream  the  first 
turbulent  spots  occur  untill  transition  is  completed.  This  scenario  is  called 
natural  transition. 

For  higher  turbulence  levels  bypass  transition  takes  place.  It  is  governed 
by  the  direct  intrusion  of  non-linear  disturbances  into  the  viscous  sublayer, 
where  they  initiate  turbulent  spots.  In  the  spot-formation  process  both  the 
length  scale  and  the  initial  strength  of  the  disturbance  may  be  important. 

Bypass  transition  is  governed  by  the  forcing  of  free-stream  disturbances 
on  the  viscous  sublayer.  To  describe  bypass  transition  the  transient  growth 
theory  has  been  developed  over  the  last  ten  years.  It  can  be  modeled  using 
the  Parabolised  Stability  Equations.  However,  this  modelling  starts  from 
the  amplification  of  waves  or  (coherent)  structures  present  in  the  boundary 
layer.  There  are  two  semi-emperical  models  ((Johnson,  1999)  and  (Mayle 
and  Schultz,  1997))  known  to  describe  the  initiation  for  bypass  transition 
starting  from  the  intrusion  of  free-stream  disturbances.  In  both  models, 
the  free-stream  disturbances  are  assumed  to  intrude  via  the  amplification 
and  subsequent  breakdown  of  laminar  modes.  The  amplification  rate  is  pre- 
scribed from  emperical  relations.  The  actual  formation  of  turbulent  spots 
(i.e.  the  onset  of  transition)  has  to  be  introduced  through  a presumed  cri- 
terion. 

These  models  result  in  the  well-known  intermittency  distributions  in- 
troduced in  (Narashima,  1957)  or  (Johnson,  1994).  In  the  derivation  of 
these  intermittency  distributions  the  influence  of  the  initation  process  only 
remains  in  the  (empirical)  parameter  xt  (the  streamwise  position  of  the 
start  of  transtion) . The  shape  of  the  distribution  is  determined  by  the  spot 
growth-  and  merge- processes.  This  means  that  the  theoretical  models  for 
natural  and  bypass  transition  both  lead  to  these  intermittency  distribu- 
tions, since  they  both  lead  to  rapid  non-linear  breakdown  and,  therefore, 
more  or  less  instantanious  spotproduction.  In  the  past  numerous  experi- 
ments have  been  performed  to  find  Rext  and  Re\  depending  on  the  pressure 
gradient,  turbulence  intensity,  compressibility,  etc. 

Recent  experiments  ((Schook  et  al.,  2001))  show  that  for  low  levels 
of  free-stream  turbulence  the  Narashima-  and  Johnson-intermittency  dis- 
tributions are  indeed  applicable.  However,  as  the  turbulence  intensity  is 
raised  to  levels  representative  for  gasturbine  flows,  the  intermittency  distri- 
bution changes  shape.  A closer  examination  of  the  experimental  conditions 
(at  different  levels  and  length  scales  of  the  free-stream  turbulence)  shows 
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that  the  Narashima-  and  Johnson-intermittency  distributions  are  not  found 
when  the  turbulence  length  scale  is  large  compared  to  the  boundary  layer 
thickness. 

The  shape  of  the  distributions,  found  at  these  length  scales,  suggests 
that  ’spots’  grow  or  shrink  depending  on  both  the  local  boundary  layer 
thickness  and  the  strength  of  the  impinging  disturbances.  This  means  that 
the  spots  must  be  created  directly  from  free-stream  perturbations.  They 
can  not  be  the  result  of  a initiation  process  within  the  boundary  layer. 
Therefore,  at  these  scales  the  mechanism  of  direct  intrusion  of  free-stream 
disturbances  dominates  the  initiation  process. 

This  raises  the  question  of  boundary  layer  selectivity;  which  properties 
of  a main  stream  disturbances  are  descisive  for  the  production  of  a turbulent 
spot.  In  the  present  study  fully  compressible  direct  numerical  simulation  is 
used  to  study  this  transition  mechanism.  In  this  paper  a numerical  study 
is  presented  in  which  a spatially  developing  laminar  boundary  layer  on  a 
flat  plate  is  subdued  to  large  scale  (non-linear)  disturbances. 

3.  Governing  equations 

For  a compressible  flow  of  an  ideal  gas,  conservation  of  mass,  momentum 
and  entropy  can  be  written  as 

0 

-Vp-V-r  (1) 

—V • q — r : Vv 

with  time  t,  density  p,  the  velocity  vector  v,  entropy  s,  pressure  p,  stress 
tensor  r,  temperature  T and  heatflux  vector  q.  The  set  of  equations  are 
closed  with  the  constitutive  relations  T = f(p,s ) and  p = pRT. 

The  use  of  the  entropy  s has  the  advantage  that  its  conservation  equation 
does  not  contain  any  acoustic  waves.  Therefore,  only  mass  and  momentum 
conservation  need  to  be  resolved  on  an  acoustic  time-scale. 

4.  Numerical  scheme 

The  equations  are  solved  with  a higher  order  finite  difference  method  on  a 
collocated  grid.  The  time  integration  uses  a third  order  low  storage  Runge 
Kutta  scheme,  which  requires  only  one  additional  storage  location.  A com- 
pact (implicit)  sixth  order  scheme  due  to  (?)  is  used  for  the  spatial  dis- 
cretisation of  all  non-advective  terms.  The  scheme  uses  five  grid  points  and 
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can  be  evaluated  with  order  N manipulations  for  N derivatives.  Advection 
is  approximated  with  a fifth  order  explicit  upwind  scheme.  In  this  way  a 
stable  integration  is  obtained  in  which  no  additional  filtering  is  required. 

Characteristic  boundary  conditions  are  implemented  according  to  the 
NSCBC  approach  of  (?).  In  this  approximation  local  one  dimensional  waves 
are  used  for  acoustic  and  convective  transport  and  supplemented  with  vis- 
cous conditions. 

Sub-stepping  is  used  in  which  acoustic,  convective  and  viscous  trans- 
port are  evaluated  in  splitted  (third  and  second-order)  time  evolution.  The 
properties  of  the  numerical  scheme  are  such  that  it  is  ideal  for  parallel 
processing.  Furthermore,  the  code  is  optimised  for  memory  access,  caching 
and  load  balancing.  The  calculations  presented  in  this  paper  typically  use 
about  400  CPU  hours  when  run  on  16  R10k3-processors  (250  Mhz). 

5.  Problem  definition 

Initially  a Blasius  profile  ( iq,(y ),  Mac*, =0.1)  is  prescribed  for  Res » equal  to 
400  (with  8 the  momentum  thickness,  this  means  that  the  Reynolds-numb er 
based  on  the  displacement  thickness  is  approximately  1000).  At  the  inflow 
boundary  the  wall  normal  velocity  v is  perturbed  (similar  as  in  (Breuer 
and  Landahl,  1990),  who  used  a spatially  introduced  streamwise  vorticity) 
according  to: 

Af(y)(  1 - 202)(1  - 2£2)  exp(— 02  - £2  + 1/2) 

- 1-~-  0 = nS*/Uoo  to  = 6 9 

z/L*  - 1/2 

< c 8 

«&(y)(i  -y/Ly)2 

The  computation  domain  has  a height  Ly  of  204*,  the  length  and  width 
are  choosen  6 Ly  and  3 Ly,  respectively.  The  lateral  length  scale  agrees  with 
simulations  by  (Breuer  and  Landahl,  1990),  who  found  that  the  result- 
ing perturbation  is  approximately  independent  of  the  spanwise  scale.  The 
typical  time-integration  continues  till  tUoo/8*=600.  In  the  simulations  a 
uniform  grid  of  193x65x192. 

The  flat  plate  is  prescribed  by  an  adiabatic  no-slip  boundary  condition. 
The  upper  boundary  is  acoustically  non-reflective.  At  the  inflow  the  velocity 
components  and  entropy  are  fixed.  At  the  outflow  the  downstream  pressure 
is  fixed.  In  the  spanwise  direction  the  domain  is  periodic. 


av 

dt 

<t>  = 

i = 
f(y)  = 
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6.  Results 

The  initial  perturbation  results  in  a checkerboard  pattern  of  (alternating) 
wall-normal  velocity.  This  leads  to  three  spanwise  rows  of  three  alternating 
negative/positive  A2-pockets.  In  the  main  stream  the  disturbance  is  rapidly 
transferred  to  the  outflow,  as  shown  in  figure  1,  leaving  behind  a pattern  of 
alternating  positive  and  negative  wall-normal  velocity,  with  a characteristic 
length  scale  of  <5*.  The  formation  of  this  pattern  goes  hand  in  hand  with 
combination  of  the  front  middle  and  middle  outer  regions  of  negative  A2  into 
an  arrow-head  structure,  shown  in  figure  2.  This  structure  is  stretched  into 
a lambda-shaped  vortex,  situated  in  the  near-wall  boundary  layer  region, 
below  36*.  The  ratio  of  leading  and  trailing  edge  velocity  of  the  lambda 
vortex  is  about  3.  Subsequently,  a rapid  local  creation  of  vortices  near  the 
legs  is  triggered,  similar  to  the  observations  of  (Singer  and  Joslin,  1994).  In 
the  later  stages  a turbulent  spot-like  structure  is  formed. 

It  is  clear  that  the  time-wise  introduction  of  the  perturbation  on  the 
wall-normal  velocity  leads  to  results  very  similar  to  those  in  (Breuer  and 
Landahl,  1990).  A difference  in  both  results  is  the  alternating  pattern  of 
wall-normal  velocity  at  tlier  center-line  which  arises  in  the  Breuer-results. 
It  is  not  clear  it  these  are  caused  by  the  difference  in  boundary  condi- 
tions ((Breuer  and  Landahl,  1990)  uses  streamwise  periodic  boundary  con- 
ditions). The  similarity  in  the  results  indicates  that  the  formulation  of  the 
initial  disturbance  (v(t)  instead  of  a spatial  stream- wise  vorticity  perturba- 
tion) is  not  important.  Instead  the  development  of  the  initial  disturbance 
into  the  A- vortex  and  the  subsequent  production  of  vorticity  at  the  stretch- 
ing legs  seems  to  be  determined  by  the  perturbation-length  scale.  Of  course, 
more  simulations  for  different  6 and  £ are  needed  to  support  this  conclusion. 

7.  Conclusions 

The  results  show  that  the  developed  code  is  suitable  for  the  study  of  the 
spot-initiation  proces.  In  the  near  future,  the  boundary  layer  sensitivity 
will  be  investigated  for  different  0 and  (.  Furthermore,  mode-selection  for 
more  random  perturbations  will  be  studied.  These  initial  disturbance  are 
placed  in  the  laminar  boundary  layer.  The  intrusion  mechanism  of  free- 
stream  disturbances  into  the  viscous  sublayer  will  be  the  subject  of  future 
studies  to  answer  the  question  of  boundary  layer  selectivity. 
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Figure  1.  Contours  of  constant  wall  normal  velocity  on  the  center  plane  at 
t = Ux/S'= 98,  245  and  392 
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Figure  2.  Top  view  of  regions  of  negative  A2  and  gray  scales  indicating  the  entropy  on 
the  flat  plate  at  tUoo/ 6*  =245,  392  and  539 


